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ABSTRACT 

We present Very Long Baseline Interferometry (VLBI) observations 217 days after the 7-ray burst 
of 2003 March 29. These observations provide further measurements of the size and position of GRB 
030329 that are used to constrain the expansion rate and proper motion of this nearby GRB. The 
expansion rate appears to be slowing down with time, favoring expansion into a constant density 
interstellar medium, rather than a circumstellar wind with an r~ 2 density profile. We also present 
late time Arecibo observations of the redshiftcd Hi and OH absorption spectra towards GRB 030329. 
No absorption (or emission) is seen allowing us to place limits on the atomic neutral hydrogen of 
Nh < 8.5 x 10 20 cm -2 , and molecular hydrogen of Nh 2 < 1-4 x 10 22 cm~ 2 . Finally, we present VLA 
limits on the radio polarization from the afterglow of <2% at late times. 
Subject headings: gamma-rays: bursts 
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1. INTRODUCTION 

Our understanding of the origin of gamma-ray bursts 
(GRBs) has continued to advance rapidly in the 
year s since the first X-rav dCosta et al.1 119971) . opti- 
cal Jvan Paradiis et al.l Il997|) and radio (|Frail et al. I 
11997(1 afterglows were discovered. In particular 
the nearby afterglow from GRB 03032 9 has s olidified 
the GRB-supernova connection l|Hiorth et alJ [2003a; 
iMatheson etafl 12003(1 , and provided th e first GRB 
with a well determined expansion rate (|Tavlor et all 
120041) . This event has pre sented a unique opportunity 
to test afterglow models (lOren. Na kar. fc Piranl [2004; 
Granot, Ramirez- Ruiz & Locb 2005), and to explore the 
environment around a GRB. 

Afterglow models invoke gas-rich environments around 
at least some of the GRBs. Several basic properties of 
this circumburst medium are presently unknown, and 
need to be addressed in order to better understand the 
afterglow and its evolution. For example, in the simplest 
emission models, the density of the ambient mediu m is 
related to the afterglow flux density l(Waxman 1997). Es- 
pecially, if GRBs exist in both gas-poor and gas-rich en- 
vironments, that could explain why afterglow is absent 
in some GRBs. Direct observations of the circumburst 
medium associated with GRBs would thus provide an 
important test of the fireball model. 

Another issue is the density profile of the medium into 
which the GRB ejecta expands. Current models pos- 
tulate a medium in which density is governed by the 
mass- loss of the progenitor star (p cx r" 2 ). At present 
however, some observations instead point to a uniform 
density in order to explain the evolution of the afterglow 
lightcurvc (Be rger et al. Il2003() . In a few cases, strong 
damped Lya absorption has been found indicating Hi 
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column densities as high as 5 x 10 21 cm~ 2 (Hiorth et al. 
I2003b() . During the expansion of the GRB ejecta, the cir- 
cumburst medium is likely to go through different stages 
of ionization. Perna & Loeb (1998) suggest that as a con- 
sequence (optical) absorption lines will vary in their EW 
with time. The rate of this variation can constrain the 
size of the absorbing region. Similarly, in the radio the 
fraction of molecular and atomic gas should be governed 
by the effective ionization, and might therefore vary in a 
similar way. 

One way to investigate the circumburst medium is via 
atomic or molecular absorption studies in the radio. An- 
other way is to measure the deceleration of the expand- 
ing fireball. In this paper we present both late time 
Very Long Baseline Interferometry (VLBI) observations 
of the afterglow of GRB 030329 and similar epoch Hi 
and OH absorption observations taken with the Arecibo 
telescope. 

Assuming a Lambda cosmology with Hq = 
71 km/s/Mpc, Q M = 0.27 and fi A = 0.73, the 
angular-diameter distance of GRB 030329 at z = 0.1685 
is d,A — 589 Mpc, and 1 milliarcsec corresponds to 2.85 
pc. 

2. OBSERVATIONS AND RESULTS 

2.1. Late Time VLBI Observations 

The VLBI observations were taken at 8.4 GHz on 2003 
November 1,217 days after the burst, with a global array 
including the Very Long Baseline Array (VLBA) of the 
NRAO 5 . Other telescopes used were the Effelsberg 100-m 
telescope 6 , the phased VLA, the Green Bank Telescope 
(GBT), the 305 m Arecibo telescope 7 , the Westerbork 

5 The National Radio Astronomy Observatory is operated by 
Associated Universities, Inc., under cooperative agreement with 
the National Science Foundation. 

6 The 100-m telescope at Effelsberg is operated by the Max- 
Planck-Institut fur Radioastronomie in Bonn. 

7 The Arecibo Observatory is part of the National Astronomy 
and Ionosphere Center, which is operated by Cornell University 
under a cooperative agreement with the National Science Founda- 
tion. 
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TABLE 1 

Observational Summary 



Date 


At 


Frequency 


Integ. Time 


BW 


Polar. 


Instrument 




(days) 


(GHz) 


(min) 


(MHz) 






20 Jun 2003 


83 


8.409 


138 


50 


2 


Y27 


1 Nov 2003 


217 


8.409 


138 


50 


2 


Y27 


1 Nov 2003 


217 


8.409 


165 


32 


2 


VLBA+EB+Y27+WB+AR+MC+NT 


30 Nov 2003 


247 


1.216 


30 


12 


2 


AR 


30 Nov 2003 


247 


1.428 


50 


12 


2 


AP, 


30 Nov 2003 


247 


1.473 


30 


12 


2 


AR 



*NOTE - EB = 100m Effelsberg telescope; Y27 = phased VLA; GBT = 105m Green Bank Telescope; WB = phased 
Westerbork array; AR = 305m Arecibo telescope; MC = Medecina 32m telescope; NT = Noto 32m telescope. 



(WSRT) tied array, and the Noto and Medecina tele- 
scopes of the Consiglio Nazionale delle Ricerche. In all 
some 56 individual and combined antennas each of 25m 
or more in diameter were employed with a combined col- 
lecting area of 0.12 km 2 . Most antennas were on-source 
for a period of 5.5 hours. All stations recorded with 
256 Mbps with 2 bit sampling in dual circular polariza- 
tion with the exception of Noto which had only a right- 
circular polarization receiver available. The observations 
were correlated at the Joint Institute for VLBI in Europe 
(JIVE). 

The nearby (1-5°) source J1051+2119 was used 
for phase-referencing with a 2:1 minute cycle on 
sourcexalibrator. The weak calibrator J1048+2115 was 
observed hourly to check on the quality of the phase ref- 
erencing. Self-calibration with a 1 hour solution interval 
was used to further refine the calibration and remove 
some slow-changing atmospheric phase errors. The fi- 
nal image has an rms noise of 24 /zJy/beam. This is 
substantially higher than the expected thermal noise of 
8/zJy/beam, in part because of partial loss of signals from 
Arecibo and Westerbork for reasons not completely un- 
derstood. 

We fit a symmetric, two-dimensional Gaussian to the 
measured visibilities on GRB 030329 and find a size of 
0.176 ± 0.08 mas. As in ITavlor eta!\ lj2004|) the error 
of the size is estimated from signal-to-noise ratios and 
from Monte-Carlo simulations of the data using identical 
(u,t>) coverage, similar noise properties, and a Gaussian 
component of known size added. The standard deviation 
of the recovered sizes, modelfitted in the same way as we 
treat the observations, was found to be 0.053 mas. 

We also obtain a position for GRB 030329 of R.A. 
10 h 44 m 49. 95955 s and Dec. 21°31T7.4377" with an un- 
certainty of 0.2 mas in each coordinate. Solving for 
proper motion using all the high frequency VLBI ob- 
servations to date, we derive fi r . a . = —0.05 ± 0.41 mas 
yr _1 and Udec. = —0.24 ± 0.41 mas yr^ 1 , or an angu- 
lar displacement over the first 217 days of 0.14±0.35 
mas (Fig. 1). Th ese observati o ns ar e consistent with 
those reported bv ITavlor et all l)2004[K and impose an 
even stronger limit on the proper motion. This limit ar- 
gu es against the cannonbal l mod el for GRBs proposed 
bv lDado. Par fc De Ruiulal pOOlj) . 

2.2. Single Dish Arecibo Observations 

Single dish observations of GRB 030329 were carried 
out with the L-band Wide receiver of the 305 m Arecibo 



Radio Telescope on 2003 November 30, for a total of 2 hr. 
The simple position-switched observations utilized all 
four interim correlator boards to observe various L-band 
frequencies. The bandwidth of each board was 12.5 MHz. 
While the first two boards were set to observe the fre- 
quency of the redshifted A21 cm H I line, recording each 
linear polarization with 2048 spectral channels, the other 
two boards were set to observe simultaneously the or- 
thogonal polarizations at the frequencies of the redshifted 
A18 cm mainline OH transitions (1665 & 1667 MHz) and 
the 1720 MHz satellite OH transition with 1024 spectral 
channels per polarization. 

Following the editing out of data suffering from radio 
frequency interference, the total on-source integration 
time for the redshifted A21 cm H I and the 1720 MHz 
satellite OH transitions was 30 min each. For the red- 
shifted A18 cm mainline OH transitions, the on-source 
integration time was 50 min. 

3. CONSTRAINTS ON THE ATOMIC AND MOLECULAR 

GAS 

3.1. Limits on Absorption 

Figures 2, 3, and 4 show pairs of total-intensity spec- 
tra of the GRB 030329 centered at the optical heliocen- 
tric velocity that corresponds to the frequency of the 
redshifted A21 cm H I line, the redshifted A18 cm OH 
mainlines, and the redshifted A18 cm OH satellite line at 
1720 MHz, respectively. 

The top spectra in Figs. 2, 3, and 4, are Hanning- 
smoothed with spectral resolutions of 3.52 km s _1 
(12.2 kHz), 5.99 km s" 1 (24.4 kHz), and 5.80 km s" 1 
(24.4 kHz), respectively. The rms noise level in each 
of these spectra is 1.02, 0.32, and 0.57 mJy beam -1 , 
respectively. The bottom plot in each figure is a five- 
channel smoothed version of the respective top spectra. 
The velocity resolution in these spectra is 8.80, 14.98, and 
14.50 km s , and the rms noise level of each spectrum 
is 0.482, 0.224, and 0.282 mJy beam" 1 , respectively. 

No emission or absorption is seen in any of the Arecibo 
spectra. We can derive a 3cr limit on the Hi opacity 
of r < 0.53 with a velocity resolution of 8.8 km s . 
This corresponds to a Hi column density limit of Nh < 
8.5 x 10 20 cm~ 2 assuming a spin temperature of 100 K, 
and uniform coverage. If the spin temperature is higher 
or the line is shallow and wide, then this limit could be 
higher. Given the extreme energies involved in the GRB, 
enough to cause a su dden ionospheric dis turbance in the 
Earth's atmosphere ijSchnoor et a/.ll2~003|) some 740 Mpc 
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Fig. 1. — The positions derived from the observations in six epochs relative to the first determination on April 1st at 8.4 GHz. Observations 
at multiple frequencies at a given epoch have been plotted separately since they are independent measurements. A circle with a radius 
of 0.2 6 mas (2a) is shown to encompass all measurements except those taken at 5 GHz, which suffer from systematic errors or et all 

2004). Taken together these observations provide a constraint on the proper motion of 0.14 ± 0.35 mas over 217 days. 



away, it could well be that nearly all the Hi along the 
line-of-sight in the host galaxy has been ionized and not 
yet recombined. Further research in this area is needed. 

In a similar manner we can place a 3cr upper limit on 
the 1667 MHz OH mainline opacity of r < 0.11 with a 
velocity resolution of 15.0 km s _1 . This corresponds to 
an OH column density limit of N OH < 1.4 x 10 15 cm~ 2 
assuming a uniform covering factor, and an excitation 
temperature of 10 K. We can further deduce a limit on 
molecular hydrogen using th e relation Nh 2 ~ 10 7 Noh 
(jKanekar fc ChengalurJ200^ . of N H . 2 < 1.4 xlO 22 cm" 2 . 

The continuum flux density of the GRB 030329 at 
1.4 GHz at the time of these observations was ~ 3.5 mJy. 
The various spectra reported here show a higher contin- 
uum flux density, because of a 15 mJy continuum source 
located at an angular distance of about 3 arcmin from 
the GRB 030329, i.e., within the primary beam of the 
Arecibo radio telescope. 

3.2. Angular Size Measurements 

Our measured size of 0.176 ± 0.08 mas is quite close t o 
the size of 0.172 ± 0.043 found bv iTavlor eta!\ ()2004|) . 
indicating a possible slowing of the burst at late times. 
This trend was already apparent from the average ex- 
pansion velocities derived from the angular size mea- 
surements on April 22 and June 20. The entire his- 
tory of expansion for GRB 030329 is shown in Fig. 5. 
The first measurement at 15 days comes from a model- 
dependent estimate of the quenching of the scintillation 
(jBerger et al. 1 12003]) . The late time curvature may in- 
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Fig. 2. — Redshifted A21 cm H I spectra toward the 
GRB 030329. The top spectrum is Hanning-smoothed with 
3.52 km s — 1 (12.2 kHz) spectral resolution. The bottom plot is a 
five-channel smoothed version of the top plot, with a velocity reso- 
lution of 8.80 km s _1 and an rms noise level of 0.482 mJy beam . 



dicate that between 83 and 217 days, GRB 030329 has 
transitioned into a non-relativistic expansion. Alterna- 
tively, it is possible (but somewhat contrived) that the 
intrinsic surface brightness profile has changed in a way 
to compensate for the expansion. For all epochs we fit 
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FlG. 3. — Redshifted A18 cm OH mainline spectra toward 
the GRB 030329. The top spectrum is Hanning-smoothed with 
5.99 km s — 1 (24.4 kHz) spectral resolution. The bottom plot is a 
five-channel smoothed version of the top plot, with a velocity reso- 
lution of 14.98 km s _1 and an rms noise level of 0.224 mjy beam -1 . 
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Fig. 4.— Redshifted OH A18 cm 1720 MHz satellite line spectra 
toward the GRB 030329. The top spectrum is Hanning-smoothed 
with 5.80 km s" 1 (24.4 kHz) spectral resolution. The bottom 
plot is a five-channel smoothed version of the top plot, with a 
velocity resolution of 14.50 km s — 1 and an rms noise level of 
0.282 mjy beam" 1 . 



a two-dimensional Gaussian model since our resolution 
does not permit us to discriminate between a Gaussian, 
ring, disk, or more complex profile. 

A gamma-ray burst drives a relativistic blast wave 
into a circumburst medium of density p whose ra- 
dius R is related to the energy of the explosion ap- 
proximately by E ~ i? 3 pc 2 7 2 , where 7 is the bulk 
Lorentz factor of the fireball. Mor e comp lete treatments 
are g i ven bv llBlandford fc McKeelll976t iChevalier fc Lil 
|2000[ iGranot fc Loebll2003|) . The density profile of the 
circumburst medium is generally taken to be either a 
1/r 2 wind, or a constant density ISM, such as one might 
find beyond the termination shock of the p rogenitor's 
stellar wind ijChevalier. Li. fc Fransson| |2*0041. The de- 



cel eration at late times currently favors the ISM model 
of l)Granot. Ramirez-Ruiz fc Loebll2005(l which shows a 
break in the expansion rate (see their Figures 4 and 6). 
Unfortunately, the large uncertainty in our measurement 
of the size of GRB 030329 at late times does not permit 
u s to discriminate strongly betwee n the various models. 

IChevalier. Li. fc Franssonl l)2004|) place the termination 
shock of the Wolf-Rayet wind at a distance of 0.4 pc 
from the progenitor. Beyond that they predict a fairly 
constant density out to the red supergiant shell at a ra- 
dius of 1.7 pc. With a current diameter of 0.5 pc for 
GRB 030329 (radius of 0.25 pc) at day 217 it is near 
the termination shock, especially if the progenitor had 
a shorter lifetime, or a relatively high density ISM has 
stalled the shock. 



3.3. Implications of low Polarization 

From 8.4 GHz VLBA observations on April 6 
iTavlor et all l)2004j) derived a 3cr limit on the linear po- 
larization of 0.16 mjy/beam, corresponding to a limit on 
the fractional polarizatio n of <1.0%. I n a co ntempora- 
neous optical observation lGreinere£ al\ (|2003(1 measure a 
polarization of 2.2+/-0.3%. The decrease in polarization 
at lower frequencies has been explained as the result of 
the source being optically thick at 8.4 GHz at these early 
times since the maximum degree of linear polarization of 
an optically thick synchrotron source is 12% while the 
maximum polarizati on of an optically thin synchrotron 
source is about 80% ljPacholczvklll970l) . To look for any 
change in polarization with time as the source transi- 
tions to optically thin at 8.4 GHz we have analyzed the 
very sensitive, late time phased VLA observations for po- 
larimetry. We find no detection at either epoch and place 
3(7 limits of <1.8% and <4.7% on 2003 Jun. 20 and 2003 
Nov. 1 respectively. The total intensity measured with 
the VLA on these epochs is 3.11 ± 0.03 and 0.75 ± 0.02 
mjy respectively. 

The emission mechanism for the GRB afterglow is 
widely accepted to be synchrotron radiation, which is 
intrinsically linearly polarized if there is an ordered com- 
ponent to the magnetic field, or if the magnetic field 
is ge nerated at the internal shocks ijMedvedev fc Loebl 
1999). The recent detection of 80 ± 20% po- 
larizatio n in the prompt 7-ra y emission from GRB 
021206 l |Coj3urri ^ Boggs 2003) , although controversial 
l|Rutledge fc Foxl 120041: iWigger et all 12004]) . has led to 
an increased interest in modeling the time and fre- 
quency de pendent behavior of the polarization from the 
afterglow (Nakar. Piran fc Waxmanll200l: lG rano il 12001 
IGranot fc Tavk>rll2005|) 

To have an observed polarization of <5% in the late 
time afterglow (see §2.2) could be expl ained as the result 
of hi ghly disordered magnetic fields ijGranot fc Konigll 
120031) . 

Propagation effects can also reduce the intrinsic linear 
polarization below detectable levels. A Faraday screen 
produced by ionized gas and magnetic fields can cause 
gradients in the observed polarization angle across the 
source, leading to depolarization if the resolution element 
of the telescope, or the size of the source if unresolved, is 
large compared to the gradients. The rotation measures 
{RM) can be related to the line-of-sight magnetic field, 
5||,by 
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Fig. 5. — The apparent expansion of GRB 030329 der ived f rom m easurements and limits on the angular size as a function of time. 
The two upper limits at 7.7 and 51 days are from lTavlor et all 12004) . as are the measurements on days 2 5 and 83 (open squa res). The 
measurement on day 15 (filled circle) is a model dependent estimate based on the quenching of scintillation by Bcrgcr et al. 1 12003ft , Finally, 
the measurement on day 217 (open square) comes from this work. 



RM = 812 



-B||dZ radians m 



(1) 



where Bu is measured in mG, n e in cm , dl in pc, and 
the upper limit of integration, L, is the distance from 
the emitting source to the end of the path through the 
Faraday screen along the line of sight. We can estimate 
the minimum magnetic field needed to produce a gra- 
dient across the source. To produce a 90° rotation at 
8.4 GHz requires a RM of 1200 rad m -2 . Assuming 
a density within the red super giant shell of 0.2 cm~ 3 
Chevalier. Li. fc Franssonll2004j) and a path length of 2 



pc, then the magnetic field strength required is B 



mG. This field strengt h is simi l ar to the equipartition 
field strengths iDouehertv et al.1 l)2003[) have found from 
modeling the radio emission from colliding-wind Wolf- 
Rayet (WR) binary systems, although the field strengths 
in these systems have probably been enhanced by the 
collision. Faraday screens in the GRB environment are 
considered in more detail bv lGranot fc Tavlorl p005|) . 

4. CONCLUSIONS 

While no detection of atomic (Hi) or molecular (OH) 
material is found towards GRB 030329, the limits of 
Nu < 8.5 x 10 20 cm~ 2 are not particularly constrain- 
ing. Observations of more shielded molecules like NH3 
towards future bursts at early times would be of interest, 
and given the expected strength of a nearby afterglow of 



~50 mjy or more, could yield detections, or place in- 
teresting limits on the amount of material in the GRB 
environment. 

Although GRB 030329 has faded considerably, it may 
still be detectable with VLBI techniques. Even crude es- 
timates of the size could differentiate between the predic- 
tions of wind and constant-density environments at these 
late times. Rad io re-brightening o f GRB 0303 29 has 
been p redicted bv lGranot fc Loebl l|2003j) . and lLi fc Sond 
(2004) estimate a level of 0.6 mjy 1.7 years after the 
burst. This re-brightening might occur as the coun- 
terjet becomes non-relativistic and therefore radiatively 
isotropic. This brightening will be accompanied by a 
temporary rapid growth in the size of the source as two, 
well-separated jets become visible, and by a change in 
the light centroid ijGranot fc Loebl l2003). If such a re- 
brightening occurs then a precise size estimate at late 
times becomes readily achievable with existing facilities. 

Rather surprisingly, we find no detectable linear po- 
larization from GRB 030329 at cm wavelengths to limits 
as low as 1%. This could indicate less order than ex- 
pected in the magnetic fields of the external shock that 
drives the afterglow, or a Faraday screen that depolar- 
izes the radio emission. All of these limits are from 8 
days or more after the burst. Given the RHESSI result 
ijCoburn fc Boggs| |2003) of large polarization from the 
7-rays, and predictions of some fireball models, it would 
be well worth searching for polarization from the prompt 
optical and radio emission (e.g., Granot & Taylor 2005). 
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In the near future the Swift satellite should dramat- 
ically increase the number of GRBs with measured red- 
shifts, revealing some that are nearby. In future VLBI 
studies of GRB afterglows at redshifts less than 0.1 it 
should be possible to image the structure of the after- 
glow. Fireball models of heating by a single relativistic 
shock front pre dict that at late times the fir eball should 
look like a ring ijGranot. Piran fc Sarilll999fl . 



see http://swift.gsfc.nasa.gov 
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